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't''  OO  abstract.  In  the  interest  of  obtaining  a  comprehensive  view 
penetration  of  bodies  into  solid  plates  some  of  the 
more  prominent  theoretical  relations  used  to  describe  pene¬ 
tration  processes  are  outlined  and  correlated.  The  velocity 
ranges  to  which  these  relations  apply  and  the  assuiaptions 
used  in  the  analytical  developments  are  summarized.  Elastic 
behavior  processes  in  impact  are  first  discussed  in  order  to 
establish  lower  bounds  on  the  impact  velocity  for  the  for¬ 
mation  of  a  crater.  Various  analytical  representations  of 
crater  fontiation  processes  are  suiuiaarized  and  the  advantages 
of  the  theoiy  of  constant  dynamic  flow  pressure  are  demon¬ 
strated  by  calculated  results  which  agree  well  with  experi¬ 
mental  penetration  data  from  several  sources .  The  nydro- 
dynarnic  theory  of  penetration  and  Opik's  theory  of  pene¬ 
tration  are  discussed  in  relation  to  these  conclusions  on 
crater  formation. 
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IMTRODUCTICJN 


The  process  of  penetration  of  particles  or  bodies,  solid  or  liquid, 
into  solid  plates  of  large  extent  is  extremely  complicated  and  is,  gener¬ 
ally  speaking,  still  relatively  little  understood.  A  rational  description 
of  the  phenomena  encountered  is  thus  difficult,  particularly  in  view  of 
the  fact  that  some  of  the  experimental  evidence  obtained  from  vlrtxially 
identical  tests  of  different  investigators  leads  to  widely  divergent 
results.  Furtheivore,  it  has  been  clearly  demonstrated  that  different 
phenomena  occ\ir  for  different  ranges  of  Inltied  impact  velocities  and 
that  different  mechanisms  of  energy  conversion  play  a  domlneuit  role  in 
these  various  domains.  Consequently,  it  cannot  be  expected  that  a  single 
comprehensive  theory  can  be  readily  evolved  which  will  predict  the  pene¬ 
tration  process  in  its  totality  for  all  materials,  iaipact  velocities,  and 
collision  geometries.  Rather,  most  Investigators  have  concentrated  on 
developing  theoretical  relations,  some  rational,  many  enqpirlcal,  deduced 
from  available  data  and  dimensional  analysis,  which  are  presumed  to  apply 
only  within  definite  limits  of  the  initial  velocity. 

It  is  the  pvuTpose  of  the  present  report  to  indicate  some  of  the  more 
prominent  theoretical  relations  determined  in  this  manner  and  to  correlate 
this  information.  An  indication  will  also  be  given  of  the  nature  and 
limits  of  the  velocity  ranges  mentioned  above  and  to  discuss  the  various 
assuiiq>tlons  inherent  in  the  analytical  development.  Some  siinple  deri¬ 
vations  are  slLso  Included  which  purport  to  apply  to  the  normal  impact  of 
spheres  and  cylinders  on  flat  plates  in  a  relatively  low  iiqpact  velocity 
regime. 


ENERGY  PROCESSES  FOR  PENETRATION 


The  energy  processes  encountered  in  penetration  over  the  entire 
spectrum  of  impact  velocities  consist  of  the  following  types: 

1.  Energy  of  elastic  wave  propagation,  E  ,  resulting  in  elastic 
(recoverable)  deformations  of  striker  and  target. 

2.  Energy  of  plastic  wave  propagation  or  plastic  flow,  Ep,  resulting 
in  permanent  deformation  of  striker  and  target. 


3.  Energy  of  comminution  or  fragmentation  of  the  colliding  bodies,  E, 
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4.  Energy  of  bulk  compression  of  the  solids,  W. 

Shock  heating  energy  due  to  the  difference  between  adiabatic  and 
shock  compression  of  the  solids,  H. 

6.  Energy  of  vaporization  emd  fusion  of  striker  and  target  material, 
Ev- 

7.  Energy  of  polymorphic  transitions  and  phase  changes  of  the 
materials, 

An  energy  balance  for  the  penetration  process  can  now  be  obtained  by 
equating  items  1  -  7  to  which  represents  the  difference  between  the 
initial  and  terminal  bulk  kinetic  energy  of  the  projectile-target  system. 
In  addition,  two  other  eqiiations  are  derived  from  the  laws  of  conser¬ 
vation  of  mass  and  momentum  of  the  system.  Provided  sufficient  Infor¬ 
mation  is  available  concerning  the  mechanical,  physical,  and  thermo¬ 
dynamic  behavior  of  the  solids  Involved,  these  equations  suffice,  at  least 
in  theory,  for  a  conqjlete  description  of  the  behavior  of  the  system. 

In  actual  practice,  however,  severe  difficulties  are  encountered  in 
the  solution  of  these  eqtiations.  In  fact,  It  is  virtually  impossible  to 
predict  damage  parameters  such  as  crater  depths  from  these  equations 
alone,  and  consequently  recourse  is  almost  universally  taken  to  the  use 
of  models  for  this  type  of  prediction.  Furthermore,  it  is  customary  to 
delete  terms  in  ranges  where  the  contribution  from  the  respective  mecha¬ 
nism  is  expected  to  be  very  small  relative  to  other  terms.  This  is 
illustrated  in  the  sequel  where  various  domains  of  the  penetration  process 
are  discussed  in  greater  detail. 


ELASTIC  COLLISIONS 


The  process  of  elastic  collisions  does  not,  by  definition,  involve 
any  permsuient  deformations  in  the  colliding  solids  and  as  such  would  not 
appear  to  be  of  basic  interest  for  the  theoretical  analysis  of  crater 
formation.  The  process  is  important,  however,  for  the  establishment  of  a 
lower  bound  of  impact  velocities  for  the  generation  of  permanent  craters 
which  can  then  sei^re  as  a  criterion  for  the  presence  or  absence  of  such 
deformations .  In  the  elastic  regime,  two  distinct  cases  arise  which  are 
treated  somewhat  differently  on  a  theoretical  basis:  (l)  the  elastic 
collision  of  two  bodies  with  plane-ended  surfaces  in  a  direction  normeil 
to  these  surfaces,  emd  (2)  the  elastic  collision  of  two  bodies  with 
cxurved  surfaces  at  the  point  of  impingement.  Both  cases  are  well-known 
(Ref.  1  and  2)  and  only  the  most  important  results  will  be  cited  here. 
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yorwifti  Collision  of  ivo  Plane-ended  Bodies 

Consider  the  collision  of  two  uniform  plane-ended  rods  of  density  f>^ 
and  p2;  cross-sectional  area  Ai  and  Aq,  and  rod  wave  velocity  c^  and  C2, 
respectively,  with  a  relative  impact  velocity  Vq,  as  shown  in  Fig.  1. 


.A,,C| 


. n 


FIO.  1.  normal  Collision  of  Two  Plane- 
ended  Bodies. 


% 


The  rod  wave  velocity  c  is  defined  here  as 

c  ./r  (1) 

where  E  is  the  modulus  of  elasticity,  and  represents  the  rate  of  propa¬ 
gation  of  a  plane  wave  in  a  relatively  narrow  rod.  At  the  instant  of 
impact,  two  elastic  waves  are  generated  in  both  bars,  whose  anplitudes 
may  be  determined  from  the  well-known  equation 

0  =  pcv  (2) 


and  the  condition  of  equality  of  forces  and  particle  velocities  v  at  the 
contact  surface,  given  by 


as 


°1  \  ~  °2  ^ 


1 

D 

°2  - 

r"2  .  ^-1 

PlCi 

_PlCi  P2C2_ 

(3) 

(4) 


The  particle  velocities  induced  in  the  bodies  are  given  by 


’a  -  .  AgPjCg 


(5) 
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For  normal  contact,  the  vave  in  the  smaller  rod  Is  conipletely  plane;  the 
vBve  In  the  larger  rod  is  Initially  spherically  divergent,  but  will 
rapidly  approach  the  character  of  a  plane  vave,  vith  an  amplitude  given 
by  Eq.  4,  provided  the  ratio  of  areas  of  the  rods  is  not  too  great.  If 
one  of  the  colliding  bodies  approaches  the  dimensions  of  a  semi-infinite 
target  (A2/A]^  »  1)>  the  vave  in  this  body  vill  be  essentially  spherically 
divergent  throughout  and  vill  travel  vith  dilatational  velocity  T2  “ 
i(\  +  20)/p  vhere  X.  and  G  are  the  Lame  constants .  The  amplitude  at  any 
distance  r  from  the  contact  point  for  this  spherical  vave  is  approximated 
by 


*^2 


1  Pl"lP2Vo 

r  p^c^  +  PgTg 


(6) 


If  body  2  cam  be  regarded  as  completely  rigid  (E2  =  cd),  the  stress 
produced  in  bar  1  can  be  determined  from  Eq.  4  vith  C2  »  00,  or  from  Eq. 
2  with  V  replaced  by  Vq;  this  sltmtlon  yields  the  maximum  stress  in  bar 
1  for  this  type  of  impact.  The  stress  at  the  contact  surface  is  main¬ 
tained  until  relieved  by  a  reflected  vave  emanating  from  a  free  end  of 
either  bar. 


The  threshhold  of  plastic  deformation  is  attained  vhen  the  stress 
produced  by  the  collision,  given  by  Eq.  4,  reaches  the  dynamic  compressive 
yield  stress  of  the  material.  This  quantity  is  usinlly  in  excess  of  the 
corresponding  static  value  by  various  eunounts;  for  mildly  strain-hardening 
materials  vithout  a  definite  yield  point  such  as  aluminum  and  some  alloy 
steels,  the  excess  is  from  about  1^  to  25%.  For  exaiiq>le,  2024-T4  aluminum 
alloy,  which  has  a  static  yield  stress  of  about  ^^,000  psi,  vovild  have  a 
dynaMc  yield  stress  of  about  65,000  psi.  On  the  other  hand,  for  materials 
with  a  definite  yield  point,  such  as  mild  steel,  the  excess  may  be  from  100 
to  200^.  Specific  literature  references  for  the  materied  under  consider¬ 
ation  should  be  consulted  for  the  determination  of  the  dynamic  coinpressive 
yield  stress  (Ref.  3-5) • 

Collision  of  Bodies  With  Rounded  Contact  Surfaces 


The  Inqplngement  of  two  bodies,  whose  contact  surfaces  ceui  be  i-epre- 
sented  analytically  by  equations  of  the  second  degree,  has  been  invariably 
analyzed  by  means  of  the  Hertz  law  of  contact.  This  relation  was  origi¬ 
nally  developed  on  the  basis  of  an  electrostatic  analogy  to  specify  the 
compression  of  two  clastic  bodies  in  contact  under  static  loeidlng  con¬ 
ditions,  but  was  edso  used  by  its  originator  emd  numerovis  subsequent 
investigators  for  a  description  of  collision  processes  involving  such 
bodies  (Ref.  1,  6,  j).  The  excellent  agreement  observed  between  the 
results  of  numerous  impact  experiments  and  the  predictions  of  the  Hertz 
theory  has  provided  confidence  in  the  postulate  that  the  Hertz  law  applies 
also  under  dynamic  conditions  provided  the  elastic  limit  of  the  material 
is  not  transcended.  Strangely  enough,  it  was  also  found  that  the  Hertz 
law  appears  to  predict  approximately  a  number  of  collision  parameters, 
such  as  the  duration  of  contact  and  the  maximum  force  of  the  collision 
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even  when  permanent  craters  are  produced  --a  condition  which  is  in  con¬ 
tradiction  with  the  basic  postulate  of  the  theory.  However,  the  elastic 
analysis  can  not,  of  course,  predict  the  dimensions  of  such  a  permanent 
indentation  since  it  presumes  completely  elastic  restitution. 

Briefly,  the  force  of  conqpression  F  between  the  bodies  is  given  by 
the  relation  (Ref.  l) 

3/2  ,  . 

F  =  (7) 


where  a  is  caLLled  the  approach  between  the  bodies,  defined  as  the  differ¬ 
ence  in  the  distance  between  two  points,  one  in  each  body,  located  along 
the  common  normal  at  the  initial  point  of  contact  and  sufficiently  far 
removed  from  the  contact  region  as  to  be  in  an  essentieLLly  unstrained 
domain  at  all  times,  at  the  initial  instant  of  contact  and  at  any  sub¬ 
sequent  instajit  while  the  bodies  are  still  compressed.  The  quantity  kg 
depends  upon  the  elastic  constants  of  the  bodies  amd  is,  furthermore,  in 
general  a  complicated  integral  function  of  the  geometry  of  the  surfaces 
in  contact.  The  contact  area  is  an  ellipse  of  semi-axes  a  emd  b,  and  the 
pressure  distribution  p  across  this  ellipse  is  given  by 


y  <  b 


(8) 


relative  to  axes  x  amd  y  extending  along  the  semi -major  and  semi -minor 
axes  a  and  b,  respectively.  Terms  a  and  b  can  also  be  evaluated  by  means 
of  integral  relations  from  the  elastic  and  geometrical  properties  of  the 
bodies.  Contact  is  terminated  whenever  the  approach  a  returns  to  zero. 


For  the  particular  case  when  the  two  bodies  (distinguished  here  by 
subscripts  1  and  2)  are  spheres  of  radius  emd  Rg,  respectively,  or  for 
the  case  of  a  sphere  and  a  plane  surface  (Rg  =  OD  ),  the  value  of  kg  assumes 
a  particularly  simple  form.  In  this  event,  with  8  =  (l  -  ii^)/rtE,  where  ^ 
is  Poisson's  ratio 


k„  = 


2  =  5«(e^  +  Qg) 


RlRg 


1/2 


for  two  spheres 


(9) 


k„  = 


3n(9^  + 


for  a  sphere  and  a  plane  surface 


The  contact  area  is  a  circle  of  radius  a,  related  to  the  approach  by  the 
equations 


2  aR^Rg 
a  = 


R. 


aR^ 


for  two  spheres 


(10) 


for  a  sphere  and  a  plane  siirface 
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respectively,  and  the  pressure  distribution,  Eq.  8,  may  be  written  as 


where  r  is  the  radial  coordinate. 


The  collision  of  two  such  bodies  can  now  be  written  in  terms  of  the 
equation  of  motion. 


°1°2 


•  « 

a  = 


-? 


(12) 


and  substitution  of  Eq.  7  then  yields  the  relation 

a  =  (13) 

where  k^^  =  (mj_  +  m2)/m]_me.  If  body  2  can  be  regarded  as  semi-infinite 
(m2/mj^  »  l)>  the  quantity  kj^  may  be  written  simply  as  k^  =  l/mj^.  Inte 
gration  of  Eq.  13  yields 


i<d2  -  v2)  .  .  I 


(IM 


Where  Vq  is  the  initial,  impact  velocity,  and  since  d  =  0  at  the  maximum 
approach  Q^,  this  term  can  be  evailxiated  from  the  expression 


a 


m 


2/5 


(15) 


The  average  pressure  across  the  contact  circle  at  the  instant  of  maximum 
approach  is 


max 

^av  2 
jta 

max 


L  -Kj 


V5  (»i  *  "s) 


nR^Rg 


(16) 


and  the  maximvim  pressure  at  this  instant,  p,n,^y,  may  be  determined  from 
Eq.  11  as 


P 


max 


av 


F 

max 

2 

na 


(17) 


located  at  the  center  of  the  circle.  The  total  duration  of  contact  may 
be  evaluated  as 
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Davies  (Ref.  7)  cosiputed  the  stress  distribution  in  a  semi -infinite 
elastic  solid  when  loaded  by  a  steel  sphere  (Fig.  2),  and  from  this  stress 
field  determined  the  values  of  the  principal  stresses  (stresses  acting  on 
an  Infinitesimal  cubic  element  of  the  solid,  so  oriented  that  no  shear 
stresses  act  on  the  faces  of  the  cube)  throughout  the  body.  This  calcu¬ 
lation  is  necessary  if  the  state  of  incipient  plastic  yielding  of  the 
semi- infinite  body  is  to  be  ascertained,  since  all  criteria  for  this  con¬ 
dition,  postulated  thus  far  for  multi-axial  stress  conditions,  involve  the 
ma^itudes  of  the  principal  stresses .  The  two  most  commonly  used  criteria, 
which  meet  the  experimenteilly  observed  condition  that  hydrostatic  pressure 
has  little,  if  any  influence  on  yielding,  are  the  so-called  Tresca  emd 
von  Mlses  conditions,  given,  respectively,  by  (Ref.  8) 


-  Oj  =  constant 

(19) 

a 

2  2  2 
(oi  -  Og)  +  (oj^  -  Oj)  +  (og  -  Oj)  =  constant 

(20) 

where  and  oj  are  the  maximum  and  minimum  (algebraic)  principal  stresses. 
The  constant  in  Eq.  19  is  usually  chosen  as  o^,  the  yield  stress  is  simple 
tension  or  compression  (which  are  assumed  to  be  equal  for  most  metals), 
while  the  constant  in  Eq.  20  is  chosen  as  2oy^.  The  two  criteria  differ 
to  some  extent  in  the  prediction  of  the  maxiiiivim  shear  stress  imatc  for 
various  stress  combinations,  i.e.,  -Xna-^  =  l/2  oy  for  Eq.  19  and  Tmnv  = 
2oy/yT  for  Eq.  20,  but  as  indicated  by  these  values,  the  differences 
are  not  too  great.  Equation  20  has  been  foxmd  to  predict  more  accurately 
the  yield  condition  for  most  ductile  metals,  but  Eq.  19  is  simpler  to  use 
6uid  has  thus  found  considerable  favor  in  analytical  investigations. 


FIG.  2.  Stress  Distribution  in  a  Semi-infinite 
Elastic  Solid  When  Loaded  By  a  Steel  Sphere. 
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The  analysis  perfomed  by  Davies  utilized  Eq.  19  for  the  prediction 
of  an  incipient  state  of  yielding.  Dy  suitable  cooiblnatlon  of  the  prin¬ 
cipal  stresses  in  the  semi- infinite  body,  he  detemlned  that  the  maxlfflum 
difference  of  the  principal  stresses  in  such  a  body  exhibiting  a  Poisson's 
ratio  of  0.286  (corresponding  to  steel)  could  be  ejqpressed  by 

“  0.63  p^  (21) 

where  is  given  by  Eq.  17*  When  Eq.  21  is  then  equated  to  the  dynamic 
yield  stress  in  slnple  tension  or  conpresslon,  it  is  possible  to  determine 
either  the  mlnlimim  Impact  velocity  required  to  produce  plastic  deformation 
for  a  material  with  a  known  value  of  the  dynamic  yield  stress  in  tension 
or  compression,  or  else  to  experimentally  determine  such  a  limiting  stress 
by  observing  the  first  noticeable  deformation  of  an  optically  flat  surface 
subjected  to  the  inqtact  of  a  sphere  dropped  from  successively  greater 
heights.  In  either  case,  the  equation  for  this  case  may  be  detezvlned 
from  Eq.  21  and  16,  with  ■■2/°^  =  Rg  =  oo,  as 


Dynamic 


=  0.63  p, 


max 


=  0.63  X 


r  ”1  1 


k 

5 


(22) 


It  sho\ild  also  be  noted  that  the  position  of  the  maximum  difference  of 
the  principal  stresses  does  not  fall  on  the  circle  of  contact,  but  rather 
a  certain  distance  vertically  below  the  center  of  this  circle,  so  that 
plastic  flow  is  initiated  in  the  interior  of  the  body.  For  the  case 
examined  by  Davies  (Ref.  7),  corresponding  to  steel,  this  distance  is 
0.46a  where  "a"  is  the  radius  of  the  circle  of  contact. 

Davies  observed  that  hard  steel  spheres  with  diameters  of  the  order 
of  0.23  inch  produce  i>ermanent  Indentations  in  armor  plate  (static  yield 
stress,  155,000  psl:  ultimate  tensile  strength,  >43,000  psi;  Brinell 
hardness  number  321)  when  dropped  from  as  low  a  height  as  0.35  cm.  This 
impact  corresponds  to  eui  initial  impact  velocity  of  Vq  =  10. 35  in/ sec  ■ 
0.861  ft/sec,  an  extremely  small  velocity,  yet  produced  a  maximum  stress 
PmaY  =  289,000  psi  at  the  center  of  the  contact  circle,  corresponding, 
from  Eq.  21,  to  a  dynamic  yield  stress  of  182,000  psi.  It  ceui  be  con¬ 
cluded  from  this  sequence  of  observations  that  If  permanent  indentations 
can  be  produced  by  solid  spheres  in  a  material  as  hard  as  armor  plate  at 
these  minimal  velocities,  that  purely  elastic  impacts  by  spheres  against 
plane  surfaces  virtually  do  not  exist  in  practice  and,  furthermore,  that 
the  elastic  range  can  normally  be  completely  neglected  under  actual  Impact 
conditions  of  this  type.  This  contrasts  sharply  with  the  situation 
encoxxntered  in  the  Impact  of  plauie  surface,  where  permanent  indentations 
are  not  expected  to  occur  wtil  the  stresses  0^  or  Og  given  by  Eq.  4 
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attain  the  dynamic  compressive  yield  stress  of  the  material,  ^Ich  requires 
substantially  larger  Impact  velocities  than  for  the  case  of  a  sphere 
c-igalnst  a  plane  surface. 


TECHNIQUES  OF  ANALYSIS  FOR  THE  FORMATIC»i  OF  PERMANENT  CRATERS 


The  two  principal  approaches  for  the  analysis  of  the  formation  of 
permanent  craters  can  be  described  on  the  following  basis:  (l)  techniques 
enqiloylng  the  theory  of  plasticity,  usually  applied  to  gui  ideally  plastic 
solid  (i.e.,  a  solid  In  which  the  yield  stress  under  uni-axlal  load 
remains  constant  with  increased  deformation;  this  contrasts  to  a  work¬ 
hardening  solid  in  which  this  stress  Increases  with  further  deformation), 
and  (2),  a  description  of  the  behavior  of  the  materials  by  meeuis  of  the 
continuity  equations  (Bernoulli  equations)  for  a  perfect  fluid,  compress¬ 
ible  or  Incompressible,  with  particular  reference  to  the  penetration  of 
targets  by  Jets.  It  should  be  noted  that  the  theory  of  plasticity  always 
postulates  complete  Incompressibility  of  the  materials.  Some  other 
theories  of  penetration  utilizing  the  concept  of  Incompressible  fluid 
flow,  such  as  advanced  by  Opik  (Ref.  9  and  10 ),  are  closely  related  to 
the  analysis  based  on  plasticity  arguments. 

The  theory  of  plasticity  has.  In  general,  been  developed  only  for 
equilibrium  processes  and  must  thus  be  adapted  to  dynamic  conditions  In 
order  to  permit  application  to  the  process  of  crater  formation.  Even  for 
static  loading,  however,  relatively  few  deformation  problems  have  been 
solved  In  a  completely  rigorous  fashion  In  the  plastic  regime.  While,  for 
example,  the  case  of  Indentation  of  a  semi -infinite  solid  with  a  plane 
surface  by  a  flat  die  of  infinite  length,  loaded  by  a  uniform  normal  pres¬ 
sure  (two-dimensional  case),  has  been  Eumlyzed  exactly  by  mews  of  the 
plasticity  criterion  given  by  Eq.  19  or  20,  (Ref.  8),  the  corresponding 
three-dimensional  problem  of  the  penetration  of  flat  circular  p\inch  into 
such  a  solid  has  been  evaluated  only  approximately  due  to  the  necessity 
of  using  a  plasticity  criterion  having  no  physical  basis  (Ref.  ll).  The 
latter,  known  as  the  Haar-Karman  criterion,  permits  an  integration  of  the 
stress  field,  which  the  criteria  expressed  by  Eq.  19  or  20  will  not  permit, 
but  requires  that  the  two  smaller  principal  stresses,  05  or  oj  are  equal 
wd  then  utilizes  Eq.  19.  The  Indentation  of  a  sphere  into  a  semi- infinite 
solid,  which  is  also  a  three-dimensional  problem,  has  likewise  only  been 
solved  by  means  of  the  Haar-Karman  criterion.  Consequently,  the  available 
solutions  of  three  -  dimens  i*  onal  problems  of  plastic  Indentations,  while 
probably  reasonably  accurate  on  the  whole,  must  still  be  regarded  as 
approximations  in  view  of  the  enforced  equality  of  two  of  the  princlpeLl 
stresses,  a  condition  which  actually  will  not  prevail. 

The  actuELl  solutions  of  all  these  plastic  indentation  problems,  regard¬ 
less  of  the  plasticity  criterion  employed,  has  been  carried  out  by  the 


9 


HAVWEPS  REPORT  7812 


method  of  slip-lines.  It  is  assvimed  here  that  plastic  flow  occurs  when¬ 
ever  the  maximum  shear  stress  reaches  a  critical  value  k  (where  k  »  l/2  ay 
according  to  Eq.  19  and  k  =  0.575  oy  according  to  Eq.  20)  and  that  the 
effects  of  elastic  strains  can  be  neglected.  The  flow  of  the  material,  can 
then  be  determined  by  integrating  along  these  slip-lines  which  are  defined 
as  the  families  of  curves  whose  directions  at  every  coordinate  in  the  body 
correspond  to  those  of  the  maximum  shear  strain  rate.  The  details  of  this 
procedure  are  quite  Involved  and  will  not  be  cited  here,  but  a  full  expo¬ 
sition  of  such  techniques  is  given  in  Ref.  11,  pp.  13^  ff. 

In  the  application  of  the  theory  of  plasticity  to  cratering  problems, 
it  is,  of  course,  necessary  to  employ  the  value  of  the  dynamic  yield  stress 
rather  than  the  corresponding  static  value  in  the  analysis.  Furthermore, 
it  is  necessary  to  account  for  the  inertia  of  the  target  material  in  the 
appropriate  kinetic  relations;  this  quantity  is  usually  neglected  in  slip¬ 
line  solutions  since  the  resistive  forces  of  the  target  normally  far 
exceed  inertial  terms  in  static  deformation  problems.  For  impact  con¬ 
ditions,  however,  the  magnitude  of  this  term  is  frequently  not  negligible. 
It  will  be  assumed  in  the  development  of  the  theory  based  upon  plasticity 
arguments,  that  fragmentation,  vaporization,  and  shock  heating  of  the 
colliding  objects  is  absent. 

From  the  available  slip-line  solution,  it  has  been  possible  to  deter¬ 
mine  the  pressxire  under  the  Indenter  which  must  be  overcome  to  produce 
plastic  indentation.  This  pressure  is  greater  than  the  yield  point  stress 
of  the  substance,  since  the  portion  of  the  semi-infinite  target  immedi¬ 
ately  below  the  indenter  is  confined  by  surrounding  material  which  must 
be  brought  to  a  state  of  incipient  plastic  flow  before  actual  penetration 
ceui  take  place.  This  is  due  to  the  fact  that  actual  chamges  in  configu¬ 
ration  of  the  target  must  be  Initiated  at  a  free  surface.  The  values  of 
the  pressura  so  computed  and  other  pertinent  information  will  now  be  given 
for  the  known  slip-line  solutions  (Ref.  6  and  ll). 

It  will  first  be  assumed  that  no  friction  exists  between  indenter  and 
target.  Then,  for  the  two-dimensional  case  of  a  flat  die  cited  above,  the 
pressure  produced  is  uniform  across  the  die  and  of  amount  p  =  2k(l  +  it/2)* 
Thus,  according  to  the  Tresca  yield  condition,  p  =  2.571  oy  and  according 
to  the  Mlses  condition,  p  =  2.96  oy.  The  slip  lines  consist  of  arcs  of 
circles  and  straight  lines  making  45-degree  angles  with  the  free  surfaces 
and  the  interface,  and  flow  at  the  interface  occurs  along  these  lines  with 
a  velocity  fy^)x  where  x  is  the  velocity  of  the  die.  This  solution  is 
exact. 

For  the  circular  punch  and  the  ball  indentation  problem,  the  solution 
is  only  approximate  due  to  the  restrictions  of  the  Haar-Karman  plasticity 
criterion.  Even  for  the  frictionless  case,  however,  the  pressure  distri¬ 
bution  across  the  interface  is  not  uniform;  the  nature  of  the  variation 
is  shown  in  Fig.  3*  However,  the  mean  value  of  the  pressure,  pgj  =  2.84  oy 
for  a  circular  punch  and  Pj,  =  2.66  oy  for  a  sphere  (Ref.  11 ),  do  not 
differ  slgniflceuitly  from  the  corresponding  value  for  the  die.  In  neither 
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case  are  the  slip  lines  uniform  straight  lines;  however,  for  the  punch, 
the  intersection  of  the  slip  lines  with  the  interface  still  occurs  at 
45  degrees,  and  thus  the  velocity  of  the  material  below  the  indenter  is 
still  given  by  (y2')x.  On  the  other  hand,  for  the  ball,  this  angle  varies 
from  about  43  to  20  degrees,  and  thus  the  velocity  of  the  material  below 
the  striker  will  range  frcxn  (y2)x  to  2.73x,  with  a  mean  value  of  about  2x. 


FIG.  5»  Pressure  Distribution  According 
to  the  Haar-Karman  Criterion  and  Slip¬ 
line  Solutions  for  (a)  a  Flat  Circular 
Punch  (b)  a  Rigid  Sphere  Indenting  a 
Semi-infinite  Target  (Ref.  ll). 


The  effect  of  friction  is  to  further  distort  the  uniformity  of  the 
pressure  distribution  across  the  face  of  the  indenter,  and  furthermore, 
to  increase  the  effective  mean  pressure  required  to  produce  plastic 
indentation.  While  a  quantitative  evaluation  of  this  effect  is  difficult, 
it  has  been  demonstrated  for  a  case  comparable  to  those  described  above, 
that  the  increase  in  pressure  is  about  (l  +  f),  where  f  is  the  coefficient 
of  friction.  This  condition  applies,  however,  only  as  long  as  fp  <  k 
where  k  =  xy  is  the  yield  stress  in  shear;  this  amounts  to  a  raaximvnn  value 
of  the  coefficient  of  friction  f  of  about  I/6,  a  very  low  value  for  the 
contact  of  metal  surfaces.  For  higher  values  of  the  coefficient  of  fric¬ 
tion,  slip  will  occur  within  the  skin  of  the  metal  below  the  indenter 
rather  than  at  the  interface,  and  the  pressure  resisting  indentation  may 
be  considerably  in  excess  of  I.5  Pra>  where  p^  is  the  flow  pressure 
opposing  penetration  under  frictionless  conditions  (Ref.  11). 
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The  above  considerations  have  been  presented  In  order  to  examine 
whether  or  not  the  concept  of  a  constant  flow  pressure  could  be  satis¬ 
factorily  employed  for  the  analysis  of  penetration  problems  based  on  a 
theory  of  plasticity  and,  If  so,  how  this  flow  pressure  could  be  adduced 
from  some  other  known  or  readily  determinable  mechanical  property  of  the 
material.  It  should  be  mentioned  that  the  mean  flow  pressure  determined 
aumlytlcally  for  the  penetration  of  a  sphere  on  the  basis  of  the  approxi¬ 
mate  Haar-Kannan  criterion  was  somewhat  dependent  on  the  depth  of  pene¬ 
tration,  decreasing  slightly  with  Increasing  depth  under  static  loading 
conditions  (Ref.  ll).  Ejqperlments  under  these  circumstances.  Indicate, 
however,  that  a  slight  Increase  of  this  parameter  with  depth  of  Inden¬ 
tation  occurs  (Ref.  ll).  Other  tests  under  dynamic  conditions  Involving 
the  Impact  of  a  sphere  on  a  plane  surface  at  relatively  low  Initial 
velocities  have  been  reported  (Ref.  12)  vrtilch  Indicate  a  virtually  con¬ 
stant  vsLLue  for  the  flow  pressure.  On  the  other  hand,  similar  experiments 
at  higher  Initial  velocities  Indicated  a  variation  of  the  magnitude  of  the 
flow  pressure  during  the  course  of  the  penetration  process  (Ref.  IJ).  It 
is  not  unlikely  that  the  flow  pressure  may  vary  with  strain  —  or  loading- 
rate.  It  should  also  be  borne  In  mind  that  the  analysis  of  the  static 
Indentation  process  was  carried  out  for  an  Ideally  plastic  material,  which 
acttially  does  not  exist,  and  that  the  effect  of  elastic  strains  was  neg¬ 
lected.  Nevertheless,  on  balance.  It  appears  as  though  the  assiunptlon  of 
a  constant  flow  pressure  does  not  violate  too  seriously  the  actual  defor¬ 
mation  pattern  of  the  material.  Furthermore,  this  type  of  postviiate  per¬ 
mits  a  drastic  simplification  of  an  otherwise  almost  intractable  problem 
and  might  even  be  rationalized  on  the  basis  of  the  current  relative  igno¬ 
rance  concerning  the  actual  material  properties  under  impact  conditions, 
which  might  lead  to  greater  errors  than  basically  Inherent  in  this 
assumption.  A  calculated  value  for  this  pressure  based  on  a  known  dynamic 
yield  stress  can  be  obtained  with  reasonable  accuracy  only  if  existing 
friction  is  virtually  eliminated;  in  this  case,  the  magnitude  of  the  pres¬ 
sure  is  approximately  given  by  ~  2.9  oy  The  presence  of  friction 
will,  however,  modify  such  a  value  significantly.  This  feature  will  be 
further  illustrated  in  a  subsequent  section  where  the  data  from  a  number 
of  experiments  employing  flat,  cylindrical  projectiles  fired  eigalnst  a 
semi- infinite  target  at  normal  incidence  are  examined. 

Any  theory  of  penetration  based  solely  on  the  concept  of  plasticity 
must  necessarily  fail  when  the  impact  velocity  attains  a  sufficient 
magnitude  to  alter  appreciably  the  density  of  either  striker  or  target. 
This  situation  probably  occurs  long  before  either  fusion  or  vaporization 
of  the  bodies  takes  place,  so  that  the  former  condition  can  serve  as  a 
limit  of  applicability  of  the  plasticity  theory.  It  may  be  possible  to 
set  a  predetermined  bound,  perhaps  of  the  order  of  1-2^  on  the  compressi¬ 
bility  of  the  substances,  beyond  which  the  use  of  a  hydrodynamic  model 
for  the  two  bodies  is  indicated.  However,  it  may  be  possible  to  use  the 
theory  of  plasticity  in  the  event  of  fragmentation  of  either  striker  or 
target  with  the  aid  of  an  appropriate  energy  balance  provided  the  other 
effects,  cited  above,  are  not  present  to  a  significant  degree. 
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Any  success  of  the  theory  of  plasticity  in  correctly  predicting  pene¬ 
tration  phenomena  is  coupled  to  the  ratio  of  the  energy  which  is  consumed 
in  the  plastic  deformation  of  the  target  as  opposed  to  the  eunount  used  in 
providing  kinetic  energy.  In  the  regime  of  low  velocity  impacts,  most  of 
the  energy  of  the  strilcer  is  utilized  in  performing  quasi-static  work 
during  the  penetration  process.  As  the  velocity  of  the  projectile  is 
increased,  a  successively  larger  percentage  of  this  energy  is  converted 
into  kinetic  energy  of  the  target  area  and  its  environs,  and  a  hydro- 
dynamic  analysis  becomes  more  realistic. 

The  hydrodyneunic  theory  of  penetration  has  found  its  greatest  use  in 
the  prediction  of  crater  depths  and  volumes  produced  in  targets  by  Jets 
(Ref.  14-19).  This  type  of  analysis  suffers  from  a  number  of  defects  on 
theoretical  grounds  alone.  The  most  serious  deficiency  is  the  assumption 
of  a  steady  stat.i  for  the  phenomenon,  a  condition  which  is  surely  not 
realized,  particularly  during  the  initial  phase  of  the  crater  formation. 
This  postulate  is  primarily  responsible  for  the  inapplicability  of  this 
theory  to  the  penetration  of  targets  by  a  single  solid  projectile,  since 
such  an  event  can  hardly  be  regarded  as  an  equilibrium  process  (unless, 
conceivably,  perforation  of  a  very  thin  target  is  considered  where  the 
change  of  velocity  of  the  projectile  is  small).  Furthermore,  the  steady- 
state  hypothesis  implies  essentially  an  infinite  reservoir  of  jet  mate¬ 
rial,  so  that  starting  and  stopping  of  the  jet  have  no  effect  on  the 
phenomenon.  This  requirement  also  implies  the  absence  of  a  velocity 
gradient  in  the  jet,  although  a  correction  for  the  existence  of  such  a 
gradient  has  been  proposed  (Ref.  15)*  A  recent  innovation  has  also  been 
eui  attempt  to  Incorporate  variations  in  jet  density  in  the  analysis. 

In  the  elementary  version  of  the  hydrodynamic  jet  penetration  process 
both  target  and  jet  are  assumed  to  be  incompressible,  which  is  an  unneces 
sary  postulate.  Vfhile  Cook  (Ref.  17)  has  added  a  correction  term  account 
ing  for  the  compressibility  of  the  target,  the  compressibility  of  the  jet 
is  not  teUien  into  account  and,  furthermore,  the  Inertial  effects  of  the 
target  are  calcxlLated  on  the  basis  of  its  original  density.  The  effect 
of  target  strength  in  resisting  penetration,  while  neglected  in  the 
origlneLL  derivation  on  the  basis  of  its  order  of  magnitude  relative  to 
the  inertial  terras,  has  been  subsequently  incorporated  in  a  nvanber  of 
derivations  (Ref.  17-19)-  The  existence  of  secondary  penetration,  i.e., 
a  further  increase  in  the  depth  of  the  crater  after  cessation  of  jetting, 
was  also  stipulated  due  to  the  necessity  of  reducing  the  inertial  pres¬ 
sure,  generated  by  the  flow  of  the  material,  to  the  level  of  the  dynamic 
yield  strength  of  the  material.  This  dissipation  process  was  postulated 
to  consist  of  radial  flow  by  Pack  and  Evans  (Ref.  l8)  whence  the  depth  of 
the  crater  due  to  secondary  penetration  was  simply  taken  as  r,  the  radius 
of  the  hole  produced  by  the  jet.  Cook  (Ref.  17)  incorporated  the  second¬ 
ary  penetration  in  the  calculation  of  the  crater  volume  by  means  of  a 
force  balance  in  which  the  secondary  penetration  process  was  governed  by 
the  saune  strength  function  as  in  primary  penetration.  The  same  author 
also  suggested  a  method  of  accounting  for  shock  heating  and  possible 
"impact  explosions"  (instantaneous  vaporization  of  the  materials).  Some 
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details  of  the  proposed  relations  covering  these  topics  will  be  repro¬ 
duced  in  a  subsequent  section. 

Another  version  of  the  normal  impact  of  an  elastic,  cylindrical  pro¬ 
jectile  on  a  plastically  deformable,  semi -infinite  target  has  been  proposed 
by  Backman  (Ref.  20  and  21 ).  Basically,  the  technique  consists  of  an 
impedance  matching  of  the  two  bodies  at  the  interface,  and  from  the  require¬ 
ments  of  equality  of  particle  velocity  and  stress  at  this  point  deriving 
relations  of  the  form  given  by  Eq.  4  and  5  for  the  stress  wave  propagation 
into  projectile  and  target.  The  return  of  the  wave  reflected  at  the  free 
end  of  the  cylinder  to  the  impact  point  produces  a  new  set  of  compatibility 
conditions  at  this  station,  which  may  be  determined  by  superposition.  It 
is  assumed  that  the  wave  velocity  in  the  projectile  is  that  of  the  rod  wave, 
c;  the  velocity  in  the  target  is  the  dilation  wave  "c,  and  wave  reflection 
in  the  target  is  ignored.  The  impact  velocity  is  restricted  to  a  regime 
which  does  not  produce  fracture  in  the  target. 

The  wave  action  is  analyzed  on  the  basis  that,  for  the  lower  range  of 
impact  velocities,  the  reaction  of  the  target  is  essentially  elastic - 
plastic  deformation  occurring  only  on  the  crater  sides.  At  higher  veloc¬ 
ities,  the  state  of  the  target  at  the  interface  is  plastic,  giving  rise 
to  a  proportionately  lower  penetration  resistance  than  in  the  elastic 
realm.  While  the  details  of  the  analysis  in  the  two  references  Just  cited 
differ  somewhat,  it  is  predicted  that  the  functional  relation  between  the 
depth  of  penetration  and  impact  velocity  consists  essentially  of  two 
straight  lines,  such  as  sketched  in  Fig.  4,  with  empirically  determined 
constants  Vq  and  vj_  selected  to  give  the  best  fit  of  the  data.  (Note; 
in  Ref.  20,  a  more  complicated  analysis  attempting  to  accovmt  for  work¬ 
hardening  of  the  solid  necessitates  the  selection  of  three  empirically 
determined  consteuits . )  The  abrupt  transition  from  elastic  to  plastic 
resistance  of  the  target  occurs  at  an  appreciable  impact  velocity, 
amounting  to  from  5OO  to  400  ft/sec  for  a  l/4- inch-diameter  projectile 
of  tool  steel  against  a  1  i/4- inch- thick  plate  of  2024-T4  aluminum. 

On  the  basis  of  the  previous  discussion,  it  does  not  appear  reason¬ 
able  that  an  elastic  reaction  should  be  expected  at  the  interface  even  at 
low  velocities,  provided  any  cratering  takes  place  at  all.  This  is  due 
to  the  fact  that  the  initial  regime  of  plastic  flow  has  been  shown  to 
occur  directly  under  the  indenter,  although  not  at  the  interface .  Wliile 
no  flow  of  the  material  is  expected  until  a  band  emanating  from  the 
initial  zone  of  plastic  deformation  has  reached  a  free  surface,  and 
while  such  a  region  of  plastic  deformation  is  confined  by  the  surrounding, 
elastically  strained  material,  it  ceui  not  be  accepted  that  the  target  just 
below  the  indenter  sho\ild  behave  elastically  in  the  domain  of  measurable 
craters,  since  slip-lines  must  always  emanate  from  this  domain  (Ref.  8). 

It  also  appears  unreasonable  that  there  should  be  such  an  abrupt  and 
distinct  transition  from  elastic  to  plastic  behavior  of  the  target  below 
the  indenter.  Rather,  it  woxild  seem  that  it  might  be  possible  to  neglect 
the  elastic  resistance  of  the  target  altogether  (except  for  the  determi¬ 
nation  of  a  minimvun  impact  velocity  producing  measurable  penetration)  and 
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analyze  the  process  on  the  basis  of  a  constant  flow  pressure.  It  will  be 
shown  in  the  next  section  that  such  a  procedure  produces  em  excellent  fit 
to  all  the  data  presented  in  Ref.  19  and  20,  as  well  as  of  additional  test 
resxxlts  provided  by  Backman.  Furthermore,  this  sin5>le  theory  requires 
only  a  single  experimentsil  constant,  the  dynamic  flow  pressure,  which  can 
be  checked  against  the  value  predicted  from  the  theory  of  plasticity. 


PENETRATION  OF  A  RIGID  INDENTER  INTO  A  SIMI- INFINITE 
TAR®T  EXHIBITING  A  CONSTANT  FLOW  PRESSURE 


It  will  be  assxjmed  for  the  purposes  of  the  present  analysis  that  the 
material  behavior  of  a  semi-infinite  target  resisting  penetration  can  be 
described  by  means  of  the  concept  of  a  constant  flow  pressvire.  It  is 
fiirther  stipulated  that  elastic  deformations  in  the  target  can  be  neg¬ 
lected  relative  to  permanent  indentations;  however,  a  threshhold  of 
permanent  deformation,  depending  upon  the  materials  considered  and  the 
initial  impact  velocity,  will  exist,  which  can  be  described  in  terms  of 
kno\m  mechanical  properties  of  the  solid.  The  striker  will  be  considered 
to  remain  entirely  elastic  throughout  the  cratering  process,  and  the 
energy  contained  in  the  waves  travelling  in  the  striker  will  be  neglected. 


First,  consider  the  normal  penetration  of  a  rigid  circular  cylinder 
into  a  semi- infinite  target  in  the  absence  of  friction  as  shown  in  Fig.  5* 
The  resistance  to  penetration  will  consist  of  the  meein  flow  pressvire,  pjj, 
and  the  inertial  resistance  of  flow  of  the  target.  The  latter  is  given 
by  l/2  pt^,  where  is  the  target  density  and  x  the  coordinate  in  the 
direction  of  penetration,  even  though  the  slip-line  flow  has  a  velocity 
of  (V7)x  at  45  degrees  to  the  interface,  since  only  the  normal  velocity 
component  can  exert  a  resistive  force  on  the  face  of  the  indenter.  The 
equation  of  motion  is  given  by 


mx 

A 


2  V 


m 


or 


X  + 


Ap^x 


AP, 


m 


2m 


m 


=  0 


(23) 


where  m  =  ALpp  is  the  mass  of  the  projectile,  and  A,  L,  and  Pp  its  cross- 

sectional  area,  length,  and  density,  respectively.  With  v  ~  x,  inte¬ 
gration  of  Eq.  23  yields 


m 

A  / 


dv 


2m 

Pt^ 


2P, 
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tan 


-1 
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2P. 


-1 


(24) 
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FIG.  4.  Initial  Impact  Velocity 

V3 .  Penetration  Depth . 


l6 


FIG.  5-  Normal  Penetration  of  a  Rigid  Circular 
Cylinder  into  a  Semi -infinite  Target  in  the 
Absence  of  Friction. 
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and 


^2p: 


dx  V  ni  . 

V  =  ^  =  /  - tan 

dt  I  p. 


m 


“"■*  >  sT  -  f  -  /  2 


../"Vt 


(25) 


A  second  integration  provides  the  result 

m  ! 
m 
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*  = 


cos 


[.  -1  sF^X.  At  w^t^n^  1  (j. 


-1  -Fx 

cos  (tan 

xn 


(26) 


emd  the  maximum  penetration,  x^^y,  occurring  at  time  t^nax  v  =  0, 
where 


,  m  .,^2  .  -1 

t  =  T  / -  ■tan  V 

max  A  i  p.p  c 

t  m 


2P, 


m 


is  then  given  by  _ 

[ P  L  /^^P+ 

X  .  J!L.i„y  i  .  _E_1„/  ^  t  1 

- .  2p_  0,  2p. 
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^t 

Vfhenever  ^  -  «  1,  Eq.  28  can  be  expressed  as 

p  '' 

p  Lv  X  D 

_  ^ _ 0  m^  _  o 


max 


2P, 


or 


m 


2p, 


ra 


(27) 


(28) 


(29) 


which  shows  that  the  dimensionless  penetration,  x^ny/L,  varies  as  the 
sqmre  of  the  initial  impact  velocity.  The  ssune  result  can,  of  course, 
be  obtained  by  neglecting  the  flow  pressure  term  in  Eq.  25  and  writing 
the  equation  of  motion  as  the  energy  balance 


1 

2 


p  Ax 
^m  max 


(30) 


It  should  be  noted  that  the  result  obtained  here  is  in  direct  contrast 

with  the  conclusions  of  Backman. 


To  account  for  the  fact  that  a  finite  initial  velocity  Vq  is  required 
to  produce  any  permanent  deformation,  the  term  Vq  in  Eq.  28  to  50  will  be 
replaced  by  the  quantity  v^  -  v^.  This  critical  velocity  v^  can  be  ccan- 
puted  from  the  elastic  shock  matching  given  by  the  modified  form  of  Eq.  4 
with  Og  replaced  by  the  yield  stress  a^,  yielding 


V 

c 


°Y  cy t  -  vpJ 

p,  c .  •  pc 

'^t  t  p  p 


(31) 
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Consequently,  Eq.  28  and  29  niay  be  expressed  as 


respectively. 

Figure  6  is  a  plot  of  approximately  100  datm  points  derived  from  the 
normal  impact  of  0.25-inch-diameter- tool-steel  cylinders  of  various  lengths 
fired  against  1  l/4-inch-thick  targets  of  2024-T4  aluminum.  The  data  for 
this  plot  were  taken  from  Ref.  20  and  21  euid  from  personal  comnunication 
from  Backman.  It  may  be  noted  that,  below  an  initielL  velocity  of 
500  ft/sec,  the  data  clearly  define  a  straight  line,  in  accordance  with 
the  second  of  Eq.  52.  In  this  regime,  the  inertia  term,  is  small  coiiq)ared 
with  the  flow  pressure,  so  that  the  approximation  inherent  in  Eq.  29  is 
valid.  Above  this  limit,  the  data  tend  to  fall  below  the  prediction  of 
the  approximation,  which  would  be  expected  in  a  logarithmically  dependent 
function.  Consequently,  Eq.  28  should  apply  in  the  domain  above  this 
velocity,  where  the  inertial  effects  become  appreciable.  In  any  event, 
the  plot  appears  to  indicate  that  the  mechanism  of  penetration  suggested 
here,  which  requires  a  velocity-squared  dependence,  or,  alternately, 
dependence  upon  the  initial  kinetic  energy  of  the  projectile,  provides  a 
superior  fit  of  the  data  relative  to  the  theory  proposed  by  Backman,  which 
indicates  that  penetration  is  more  nearly  a  function  of  the  initial  pro¬ 
jectile  momentum. 

There  remains  only  the  evaluation  of  the  constant  flow  pressure  and 
the  value  of  the  term  Vq  from  the  data  and  a  comparison  with  predicted 
magnitudes  on  the  basis  of  the  theory  of  plasticity  and  Eq.  31.  The 
values  determined  from  Fig.  6  yield  p^  =  225,000  psi  and  Vg  =  65  ft/sec. 

The  static  yield  stress  of  2024-T4  alumin\an  is  approximately  50,000  psi, 
and  the  corresponding  dynamic  value,  using  a  25^6  increase,  is  about 
62,500  psi.  Accordingly,  the  constant  flow  pressure  should  be  2.9oy,  or 
184,000  psi.  The  discrepancy  between  measured  and  predicted  value  is 
surprisingly  small,  especially  when  it  is  considered  that  the  analysis 
neglected  frictionsLL  effects  which  would  raise  the  value  of  the  predicted 
flow  pressure  considerably.  The  value  of  v^.  computed  from  Eq.  51  is 
66  ft/sec  using  the  static  value  for  oy  or  76  ft/sec  with  the  dynaunic 
value  for  oy  which  may  be  considered  to  be  excellent  agreement  with  the 
measured  value.  The  latter  is  quite  difficult  to  determine  fiom  the 
graph  and,  furthermore,  might  yield  a  low  value  of  v^  if  the  initiail 
impact  was  even  slightly  removed  from  normal.  Considering  the  approxi¬ 
mations  in  the  analysis,  and  the  uncertainties  of  material  behavior, 
agreement  between  theory  and  data  is  eminently  satisfactory,  particularly 
since  the  material  constants  ceui  be  predetermined  by  simply  measuring  the 
dynamic  yield  stress  of  the  target. 

Consider  next  the  indentation  into  a  semi-infinite  target  of  an 
undeformable  sphere  of  radius  R  and  mass  ra,  striking  with  initial  velocity 
v^,  as  shown  in  Fig.  7>  For  the  range  of  lnq)act  velocities  where  total 
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ionersion  does  not  occur,  the  equation  of  motion  In  the  absence  of 
friction  Is 

-nx  .  jca^  |pb  +  I  Pt^^]  <53) 

vhere  a  Is  the  radius  of  the  Indentation.  Prom  the  relation 

a^  -  2Rx  -  (54) 

Eq.  33  omy  be  written  as 

.2 

X  +  (2Rx  -  x^)  +  -~  *(2Rx  -  x^)  -  0  (35) 

2  2 

Making  the  substitution  s  x  »  -x5  s  -x5t"  in  Eq.  35  yields^ 

dt  dx 


t"  -  |~t'(2Rx  -  x^)  -  ^  (f)^  (2Rx  -  x^)  =  0  (36) 

1  Ay  ^ 

and  further  substitution  of  2  »  =  (^)  leads  to  the  equation 

g  +  B  (2Rx  -  x^)  +  (2lbc  -  x^)  »  0  <57) 

This  Is  a  standard  linear  differential  equation  whose  Integration  Is 
given  by 
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(2RX-X  )dx 

(38) 


From  this  equation,  the  maximum  penetration  x  with  v  ^  0  Is  determined 
from 


p  ^  (Rx^ 

1  t  o  ^  m  max 

"  2 


max 


(39) 


which  is  obviously  only  valid  for 


P..V 
t  o 

2Pn. 

01 


<  1. 


Suggested  by  Dr.  C.  J.  Thorne  of  the  U.  S.  Naval  Ordnance  Test 
Station,  China  Lake,  California. 
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FIG.  7*  Indentation  Into  a  Semi-infinite  Target 
of  an  Undeformable  Sphere  of  Radius  R  and  Mass  m. 
Striking  With  Initial  Velocity  v^. 


For  shallow  ionersion, 


«  2Rx,  Eq.  39  reduces 


to 


max 


2P, 
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(to) 


while  for  low  velocities,  when  the  inertial  term  can  he  neglected,  Eq.  33 
has  the  solution 
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In  the  event  that  the  last  term  on  the  right-hand  side  of  Eq.  3I  can  also 

be  neglected,  corresponding  to  tt  «  2Rx,  the  maximum  penetration  may  be 
approximated  by  the  relation 
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max 


=  y 
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2nRp 

ra 


V 
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(42) 


It  may  be  noted  that  the  maximum  penetration  is  here  predicted  to  vary 
linearly  with  impact  velocity,  in  contrast  with  the  predictions  for  the 
cylindriccLL  indenter  as  given  by  Eq.  32.  Such  approximately  linear 
variation  has  been  observed  by  a  number  of  investigators  for  the  partial 
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Immersion  of  hard  spheres  into  softer  metals  (Ref.  and  22).  Figure  8 
shows  some  experimental  results  from  Ref.  13  regarding  the  penetration 
of  a  0. 50- inch- diameter  hard  steel  sphere  into  the  plane  end  of  a 
0.50- inch-diameter  bar  of  2024-T4  aluminum.  According  to  this  data,  the 
mean  flow  pressure  is  calcvLLated  as  Pj,  =  325,000.  While  this  value  is 
considerably  at  variance  with  that  obtained  from  Fig.  6,  the  calculation 
is  much  more  approximate,  aind  both  friction  and  lateral  motion  of  the 
material,  including  the  plling-up  of  the  aluminum  around  the  crater  edge, 
have  been  neglected. 


0  100  200  300 

Vo  .  INITIAL  VELOCITY,  ft/sec 


FIG.  8.  Penetration  of  a  Steel  Sphere 
Into  a  2024-T4  Aluminum  Bar. 


Once  full  immersion  has  occurred,  Eq.  33  to  42  are  no  longer  appli¬ 
cable.  The  einalysis  may  then  be  carried  out  with  the  aid  of  Eq.  23  to 
32,  with  the  area  term  A  replaced  by  and  Lpp  replaced  by  m/itl^. 
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HYDRODYNAMIC  THEORIES  OF  PENETRATION 


The  slniplest  hydrodynamic  model  of  penetration  assumes  the  existence 
of  a  continuous  Jet  of  a  material  of  density  pj  and  constant  Jet  velocity 
Vj  penetrating  a  semi-infinite  target  of  density  pt  with  penetration 
velocity  U,  also  assxuned  constant.  A  pictorial  representation  for  this 
model  is  shown  in  Fig.  9  for  both  a  stationary  and  a  moving  frame  of 
reference.  For  the  assumed  steady-state  conditions,  Bemovilli's  equation 
at  the  stagnation  point  yields  the  expression  (Ref.  9>  10,  l4,  and  13) 

P,  -  i  p/  -  I  Pj(Vj  -  <lf  (U3) 

where  pg  is  interpreted  as  the  stagnation  pressure.  Equation  ^3  assumes 
that  the  Jet  pressure  is  much  greater  than  the  yield  strength  oy  of  the 
target;  more  generally,  Eq.  43  may  be  expressed  as  (Ref.  17) 

I  p^U^  +  k  =  I  Pj(Vj  -  U)^  (44) 


where,  in  the  lower  velocity  region,  k  may  be  considered  as  the  dynamic 
yield  strength  oyj)  of  the  material,  and  in  the  higher  velocity  regimes 
may  also  Incorporate  other  dissipative  effects  such  as  compression  and 
heating  of  the  target,  radiation  of  shock  waves  in  the  target,  and  energy 
losses  due  to  laqtact  explosions,  when  present.  For  a  continuous  Jet  of 
constant  length  L  under  the  assumption  that  steady  state  is  attained 
Instantaneously,  the  penetration  P  cam  be  expressed  as 


P  = 


Udt  =  Ut.  =  UL/(v,  -  U) 


m 


where  ty  is  the  penetration  time.  For  Eq.  43,  the  penetration  may  then 
be  written  as 


P  =  L 


while,  for  Eq.  44,  this  expression  becomes 


P  *  L 


1 


"J"J 


2 


Alternatively,  am  equation  accounting  for  the  yield  strength  of  the 
target  has  been  proposed  in  the  form  (Ref.  l6) 


P 


-  Ly-i  (1 


(W) 
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\A)Bre  q  is  a  function  of  pj  and  p^;  estimates  Indicate  that  the  terms 
q(or^p/p4Vj2)  may  attain  a  value  as  high  as  0.3  for  armor  plate  subjected 
to  s&ped*'charge  action,  but  Is  negligible  for  a  leaul  target. 


JET  VELOCITY  Vj 
JET  DENSITY 


777TT 


'////// 

TARGET 
DENSITY  A 

^  PENETRATION 
TARGET  VELOCITY  U 

STRENGTH  k 


FIG.  9*  Hydrodynamic  Model  of 
Penetration  for  Both  a  Station¬ 
ary  and  a  Moving  Frame  of 
Reference . 


For  a  Jet  consisting  of  particles  rather  than  a  continuous  stream, 
it  may  be  assximed  that  inelastic  impacts  are  produced  against  the  target 

surface.  The  pressure  exerted  by  the  Jet,  considered  as  the  change  of 
momentum  divided  by  the  Jet  area,  A,  Is  then  given  by 

P  =  P/(Vj  -  U)^/A  =  I  (»»9) 
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and  Eq.  46  becooes 


(50) 


It  may  thus  be  concluded  (Ref.  l4)  that  for  any  Jet  striking  a  target, 
the  expression  for  the  penetration  In  the  simple  model  described  above 
can  be  written  as 


P 


where 


1  <  \  <  2 


(51) 


when  the  strength  of  the  target  Is  neglected.  Experimental  verification 
of  the  validity  of  Eq.  51  (Ref.  19)  did  not,  however,  provide  any  indi¬ 
cation  that  the  parameter  k  ever  significantly  exceeded  a  value  of  unity. 
The  possibility  was  suggested  that  the  actual  Jet  did  not  fracture  Into 
sufficiently  small  enough  components  so  that  the  mecheuilsm  described  by 
Eq.  49  was  a  valid  description  of  the  process.  For  single  pellets  of 
comparable  size,  on  the  other  hand,  the  inelastic  linpact  process  stipu¬ 
lated  by  this  relation  has  been  substantiated  by  observation. 

The  simple  model  described  above  is  not  always  satisfactory  in  com¬ 
parison  with  experiments.  At  very  high  Jet  velocities,  the  target 
strength  can  be  neglected,  but  other  dissipative  effects  must  be  con¬ 
sidered  in  terms  of  the  quantity  k;  at  low  Jet  velocities,  the  dynamic 
strength  of  the  target  cam  not  be  disregarded.  Furthermore,  Eq.  46,  47, 
48,  50,  auid  51  do  not  indicate  any  dependence  of  the  penetration  upon 
standoff,  whereas  experimental  results  (Ref.  18)  show  a  variation  with 
this  parameter;  the  effect  of  standoff  is  primarily  the  creation  of  a 
velocity  gradient  in  the  Jet  which  bears  upon  the  fundamental  steady- 
state  assumption  in  the  derivation  of  the  relations.  It  has  been  shown 
that,  for  a  linear  velocity  gradient,  the  formulas  given  above  can  be 
applied  to  most  metals  at  a  given  stemdoff,  provided  the  latter  is  not 
too  small  (Ref.  16) . 

In  addition  to  the  considerations  cited  above,  there  will  exist  a 
further  expansion  of  the  target  material  beyond  the  instant  of  termination 

of  Jet  action  due  to  the  presence  of  a  dynaxolc  flow  pressure  causing  flow 
of  the  target.  This  additional  penetration,  called  residual  or  secondary 
penetration,  is  concluded  only  when  this  pressure  has  been  reduced  to  the 
value  of  (tjtj.  One  version  of  this  secondary  penetration  Pg,  based  upon 
the  model  shown  in  Fig.  9^  an  assumed  hemispherical  expansion  during 
residual  penetration  and  a  xmifom  pressure  distribution  at  the  contact 
sxxrface  during  primary  penetration  (Ref.  I8),  simply  equates  P2  to  the 
radius  of  the  hole  produced  by  the  Jet.  When  the  Jet  is  considered  as  a 
rod  penetrating  a  half-space,  the  plasticeilly  deformed  hole  In  the  target 
is  essentially  cylindrical  (Ref.  IT),  and  for  large  values  of  L/A  and  P/A, 
the  final  crater  volume  Vf  can  be  expressed  as 
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Vj  =  I  kdAdP/(i  p^U^)  (52) 

where  Af  is  the  final  cross-sectional  area  of  the  hole.  Upon  substitution 
from  Eq.  44  and  47  in  Eq.  32,  it  may  also  be  shown  that  the  final  crater 
volume  can  be  written  as 

Vj  =  I  LAPjVj^/4k  =  T^Ak  (53) 


where  k  is  the  average  dynamic  yield  strength  during  the  secondary  pene¬ 
tration  process  and  Tq  is  the  kinetic  energy  of  the  Jet  Just  before 
striking  the  target.  The  threshhold  of  Jet  velocity  Just  producing 
plastic  deformation  of  the  target,  virp,  may  be  approximated  from  Eq.  53 
by  assiuolng  no  flow  and  eqriating  Vf  to  the  projectile  voltane  LA,  so  that 


(54) 


Presumably,  the  models  considered  above  are  valid  in  a  region  of  Jet 
velocity  Vj  ranging  from  v;j>  to  Vq,  where  v^  is  the  criticsLL  velocity 
beyond  which  impact  explosions  may  be  expected. 

Observed  values  of  the  ratio  U/vj  as  a  function  of  vj  from  tests 
involving  the  penetration  of  shaped  charges  with  steel  liners  into  steel, 
aluminum,  and  lead  targets  (Ref.  17  and  19)  have  been  compared  with  the 
predictions  of  Eq.  4)  in  the  form 

U/v.  =  — ■  ■  -  =  constsuit  for  a  particular  set  of  materials  (55) 

The  variation  in  Jet  velocity  was  obtained  by  increasing  the  stemdoff  of 
the  charge.  It  was  found  that  Eq.  55  provides  a  surprisingly  good  approxi¬ 
mation  with  the  experiiaental  data  in  the  range  4vm  <  Vj  <  v^,  but  indicated 
a  considerable  deviation  in  the  regime  Vip  <  vj  <  4vip.  Use  of  Eq.  44  with 
an  evaluation  of  the  term  k  by  theoretical  considerations  provided  a  much 
better  correlation  in  this  lower  velocity  regime  (Ref.  17)*  Conq)arison  of 
test  results  of  the  final  crater  volume  with  Eq.  55>  using  the  static 
yield  strength  of  li,  proved  to  be  inconclusive. 

The  work  of  compression  per  unit  volume  on  the  target  may  be  deter¬ 
mined  from  the  relation 


k  = 
P 


pdV/V 


(56) 
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which  may  be  evaluated  In  conjunction  with  any  of  the  equations  of  state 
of  the  target  material  suggested  in  the  literature,  such  as  that  given 
in  Ref.  17  or  the  Pack,  Evans,  and  James  relation  (Ref.  2^) 


P 


B, 


(Pq/p) 


BgU 


”0  ) 

i 


(57) 


(with  pressiires  in  megabars  and  densities  in  g/cm^),  where  constants 
and  Bg  for  some  metals  are  given  in  Table  1  (Ref.  2k). 


TABLE  1.  Constemts  for  Eq.  ^7  for  Severed 
Metals  as  Determined  From  Hydrostatic 
Coopressiblllty  Measurements 


Material 

,  Megabeirs 

1 

! 

pQ  ,  g/cm^ 

X* 

Aluminum. . . 

0.2209 

9.802 

2.699 

Copper . 

0.3810 

10.029 

8.88 

Iron . 

0.624 

8.099 

7.88 

Lead . 

0.1002 

12.579 

11.54 

Magnesium. . 

0.5996 

5.115 

1.741 

Uranium. . . . 

0.1172 

25.008 

10.70 

It  should  be  noted  that  the  pressure  level  determining  the  contribution 
of  the  work  of  compression  as  well  as  that  for  shock  heating  should  be 
taken  as  the  flow  pressxire  of  the  target,  l/2  energy  due  to 

shock  heating,  representing  the  difference  between  shock  and  adiabatic 
compression  T(S^,  -  8,^),  where  T  is  temperature  and  S  is  entropy,  has  been 
shown  to  be  approxliaately  given  on  a  unit  volxane  basis  by  the  relation 
(Ref.  24) 


t  Mf-  -  1)^/  (f-)’  (58) 

Here  kj^  is  the  pressure  associated  with  this  phenomenon,  Cg  is  the 
average  heat  capacity  per  unit  volume,  and  Z  is  a  constant  conqiuted  as 
1.88-IO^,  1.55'10^>  1.84'10^,  and  1.26’10**  °K  for  iron,  lead,  copper, 
and  aluminvmi,  respectively  (Ref.  17).  The  value  of  k  is  thus  given  by 

'59) 
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when  the  energy  associated  vith  shock  vave  radiation  is  neglected.  A 
good  analytical  description  for  this  component  has  not  as  yet  been 
devised . 

It  has  been  suggested  (Ref.  16}  that  the  threshhold  for  lapact 
explosion  will  be  reached  vhen  the  work  of  compression  and  shock  heating 
reaches  a  critical  level  defined  by 

*  “Si'c  ■  oA/"  (*» 

where  Is  the  molar  cohesive  energy  and  N  the  atomic  weight  of  the 
material.  In  Eq.  60,  other  contributions  to  the  quantity  k  have  been 
neglected  on  an  order  of  magnitude  basis.  The  total  volume  of  the  crater 
generated  under  these  conditions  consists  of  two  zones:  (l)  a  region  of 
the  target  where  vaporization  occurs,  denoted  by  V|^,  and  (2)  the  region 
of  secondary  plastic  deformation  V^.  In  the  first  zone,  it  may  be  assumed 
that  the  conditions  are  the  same  as  if  an  explosive  having  the  same  avail¬ 
able  energy  were  detonated  in  the  same  volume  The  energy  is  then 
given  by 

’'.(“p  *  >S.>c  -  f  Vt"! 

from  which 

'V  *  5  *  ‘S.U  *  f 

2 

where  m  is  the  total  projectile  mass  and  l/2  mvj  is  its  kinetic  energy 
before  impact.  The  crater  developed  during  the  second  zone  is  again 
described  by  Eq.  so  that  the  total  volume  for  this  type  of  Impact 
ceua  be  represented  by 

''.(1  *  (65) 

Since  p+u|  »  6k,  the  crater  produced  by  raporizetion  of  the  target 
material  Is  very  small  compared  to  that  produced  by  residual  plastic 
deformation . 


OPIK’S  THEORY  OF  CRATERINC^ 


A  theory  of  penetration  of  meteorites  has  been  developed  based  upon 
a  model  considering  both  the  meteorite  and  the  tenrget  as  incompressible 
fluids.  For  convenience,  the  striker  is  assumed  to  consist  of  a  cylinder 


2 
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of  radius  Tq,  height  ho  =  2ro,  and  mass  m  =  2«Pttr^,  with  as  the 
meteorite  density.  In  actual  practice,  the  exact  shape  of  the  meteorite 
shovild  not  exert  a  deciding  influence  on  the  phenomenon.  The  penetration 
process  is  assumed  to  be  that  shown  in  Fig.  10,  with  the  striker  flattening 
at  right  angles  to  the  direction  of  motion,  while  the  target  strength  is 
characterized  by  a  pressure  p])  representing  the  resistance  to  plastic  flow 
or  fracture;  this  parameter  is  entirely  cumlogous  to  the  dynamic  flow  pres¬ 
sure  assumed  in  previous  models.  If,  now,  r  and  h  represent  the  dimensions 
of  the  striker  during  penetration,  conservation  of  mass  yields 

r^h  =  2r^  (64) 

Denoting  the  surface  coordinate  of  the  meteor  in  contact  with  the  target 
by  z  and  the  coordinate  of  the  mass  center  of  the  meteor  by  z,  then 

•  • 
h  =  2(z  -  z)  and  from  Eq.  64  r  =  ^(r/r^)^(z  -  z)  (65) 


The  equation  of  motion  for  the  meteorite, 
is 

-  2 
mz  =  -«r  p 


striking  with  initial  velocity 

(66) 


where  p,  the  penetration  resistance,  is  assumed  to  be  given  by  the 
relation 

P  =  I  +  Pp  (6t) 


and  the  reslsteuice  to  radial  expansion,  p',  is  hypothesized  as 


.  1  *2  ^ 

P  =  2  ^t^  ^  Pd 

It  is  now  assumed  artificially  that 

.  1  *2 
P  =  P-  2  V 


(68) 


(69) 


This  corresponds  to  a  Bemoxilll  equation  for  irrotational  flow  in  a 

reference  frame  moving  with  velocity  z.  Here  p,  the  Beniovilli  constant, 
would  be  the  stagnation  pressure;  but,  of  course,  the  flow  is  not 
actually  steady. 

Equations  6^4-69,  together  with  initial  conditions  r  =  rg,  z  =  0, 

T  =  Vq,  ceui  be  integrated  numerically,  as  detailed  in  Ref.  9  and  10. 

The  volume  of  the  crater  generated  can  also  be  deduced  as  a  rough  esti¬ 
mate.  It  may  be  observed  that,  in  spite  of  the  relatively  crude  ass\iii5)- 
tions  underlying  the  present  analysis,  there  exists  good  agreement  between 
the  predictions  concerning  the  crater  volume  8u:>d  measurements  of  crater 
dimensions  for  steel  jets  penetrating  steel  targets  (Ref.  lO). 
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TARGET  DENSITY 
STRENGTH  PARAMETER 


INITIAL  CONDITIONS 


FLOW  PATTERN  DURING  PENETRATION 


Po 


FIG.  10.  Flow  Pattern  During  Penetration.  Assuned 
penetration  process  for  meteorltic  iopact. 
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